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An MC@NLO for heavy quarks in photoproduction is presented. This is the first 
lepton-hadron process to be included into MC@NLO. To construct an MC@NLO pro- 
cess dependent so called MC-subtraction terms need to be calculated. The resulting 
calculation is compared to a fixed order NLO calculation and the HERWIG event gen- 
erator and is shown to perform well. The slides can be found in [1]. 

Schematic Overview of the Calculation 

MC@NLO [H [3| is a generator for calculating next-to-leading order matrix elements, com- 
bined with a Monte Carlo event generator. When constructing an MC@NLO one needs to 
match the NLO matrix clement with the parton shower of the generator. If wrongfully done, 
some of the configurations produced by the real gluon emissions in the matrix clement may 
also be produced by the first branching in the parton shower. Therefore, to avoid to dou- 
bly count these configurations the first branching in the parton shower is subtracted from 
the real emission part of the matrix element. This method is called 'modified subtraction', 
since it is based on the subtraction method of performing NLO calculations. The terms 
being subtracted are called 'Monte Carlo subtraction terms'. These terms are dependent 
upon which process is being calculated as well as on the parton shower itself. MC@NLO 
has previously been constructed for many processes, all in hadron-hadron collisions. In the 
following will be presented an MC@NLO for heavy quarks in photo production which is the 
first MC@NLO to be constructed for lepton-hadron collisions. The NLO calculation is taken 
from FMNR [3] and the parton shower from the HERWIG event generator [S] . 

For the NLO matrix element in this process, there are three kinds of contributions, 
namely the Born term, the real emissions, and the virtual corrections. The total NLO rate 
is schematically calculated as 

^NLO = J d03 (da«'=-') + J d(^2 (da^-" + daV"t"-i) ^ 

where the integral is over the phase-space where n indicates how many particles there 
are in the final state. These two integrals arc separately divergent, even though their sum is 




Figure 1: Diagrams of the first possible branching in the HERWIG parton shower. Charge 
conjugated diagrams are implied. 
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Figure 2: The pt and 77 spectra of the b-quark as given by FMNR, HERWIG and MC@NLO. 
Parameters are as in table [1] and parton densities are MRST2002NLO 



finite. To solve the integrals a term precisely canceling the divergencies is subtracted from 
the first integral and added to the second integral: 

CTNLO = J d03 (da«<='^' - da^"^^*' -) + J d(t>2 {da^°"' + da^'"'""^^ + da^"*^*' -) . (2) 

This procedure is called the 'subtraction method', ft should be noted that when integrated 
over respective phase-space the subtraction terms cancel exactly. 

In order to match such a calculation with the parton shower without doubly counting 
any configurations the subtraction method in eq. Q is modified. Each branching in the 
parton shower contributes by an order of ag This means that only the first branching in the 
parton shower could give rise to double counting, since all other branchings give a prediction 
at orders higher than NLO. The terms canceling the divergencies present in the integrals 
are then replaced by so called 'Monte Carlo subtractio terms', which also subtract the first 
branching in the parton shower from the real emissions, and the non-branching probability 
to the second integral: 

iNLO = I d<^3 (da«^-' - dag^C) + J dcj,, (da^"" + da^-'-i + daf"") . (3) 

The MC subtraction terms thus need to cancel all the intermediate divergencies of the 
calculation, and also ensure that the first branching in the parton shower is not counted 
twice, in such a way that the total NLO rate is still preserved. After integration, these term 
do not cancel exactly, as was the case in eq. ([2]), but this is compensated by the parton 
shower. 

As explained above, to construct an MCQNLO, the MC-subtraction terms have to be 
calculated. For heavy quarks in point-like photoproduction there are two types of processes 
to consider, namely: 

75 QQg 

iq/q QQq/q (4) 

The branchings neaded to be subtracted for these processes are shown in FigH] MC- 
subtraction terms for all these branchings have been calculated. 
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Process 


FMNR 


MC@NLO 


Charm [mb] 


0.8423 ±0.003 


0.8423 ±0.003 


Beauty [pb] 


5215 ± 10 


5199 ±30 



The MC-subtraction 
terms need to cancel the 
divergencies present in 
the NLO-calculation. In 
order to test this, a reg- 
ularized version of the 
MC-subtraction terms are 
made to approach the di- 
vergent soft and coUinear 
limits, and are compared 
to a regularized version of 

the matrix element in the limit. This test has been performed succesfully. 

To be sure that the double counting is prop- 
erly canceled, by definition the total rate given 
by MC@NLO must coinsidc with that from the 
NLO calculation. The predictions for charm and 
beauty production is compared in table [1] 



Table 1: Comparisons of total prediction in FMNR and 
MC@NLO for charm and beauty production. The errors are 
statistical. Here nic = 1.5 GeV, nib — 4.95 GeV and ^b. — f^F = 
''TT'c/b,t- Parton densities are cteq5m 




Figure 3: The logarithm of the combined 
Pt of the heavy quark pair as given by 
FMNR, HERWIG and MC@NLO. Also 
shown is MC@NLO without parton show- 
ers. Parameters are as in Figl^l 



Results 

The results from the new MCQNLO for heavy 
quarks in eP-photoproduction are compared 
with predictions from HERWIG and FMNR. 
Only results on parton level arc shown, and only 
from the point-like part of the calculation. The 
HERWIG predictions have been multiplied by a 
K- factor of 1.4. 

In Fig|2]the pt and rj of the fe-quark is shown. 
For these observables no large differences be- 
tween the calculations are expected, since they 
are not very dependent upon NLO effects. How- 
ever, as seen in Figl2]the MC@NLO prediction 
is always very close to that of FMNR. Small dif- 
ferences are observed between HERWIG and the 

other calculations, especially in the very forward and backward parts of the 77 distributions 
and for large pt > 10 GeV. 

An observables where NLO effects are expected to be large are the combined pt of the 
heavy quarks. By momentum conservation, this is equal to the combined pt of all other 
final state partons. The results of this observable is shown in FigO In the case of FMNR, 
for small pt this observable show the pt distribution of the light parton emitted by the real 
emission, which diverges for small pf. At pt(bb) = (not shown in the figure) the virtual 
terms instead dominate giving a negative cross-section. It is therefore evident that the 
NLO-calculation is not sufficient to describe this quantity. HERWIG, on the other hand, 
has resummation to all orders via the parton shower. But the parton shower is still unable 
to describe the hardest part of the pt-tail for the heavy quark pair. MC@NLO also has 
resummation to all orders via the parton shower, but the NLO-matrix element also manages 
to reproduce the hard tail. Shown in FigO is also the prediction of MC@NLO without 
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parton showers. This shows how the MC-subtraction is different from the NLO subtraction 
in FMNR. It should be noted that the total cross-section in all four calculations are equal. 
Another observable where NLO effects are ex- 
pected to be large is the difference in azimuthal 
angle (j) between one of the heavy quarks and the 
hardest other parton. The result for this observ- 
able is shown in Fig. [D At LO, when no light 
parton is emitted, the heavy quarks are always 
back-to-back. In FMNR, for the light parton to 
be harder than the heavy anti-quark the angle 
needs to be larger than 90°. In HERWIG there 
are almost always more than three partons avail- 
able and thus may the hardest other parton have 
any angle relative to the heavy quark. As seen 
in the figure, this is also the case for MC@NLO. 




A<Kb, hard) ["] 



Conclusions 



Figure 4: The difference in azimuthal an- 
gle (j) between the 6-quark and the hardest 
other parton. 



The first version of MC@NLO for heavy quarks 
in photoproduction has been presented. The 
MC-subtraction terms have been constructed 
and successfully tested. Also, the predictions at 

parton level agree with that of FMNR whenever NLO is expected to be important, and 
with that of HERWIG whenever the parton shower is expected to be important. The full 
MC@NLO for heavy quarks in photoproduction is expected to be presented soon, including 
parton showers and hadronization. This will be a tool for comparison with data from the 
HERA experiments. 
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